Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

NDUFS7 variant in dogs with Leigh syndrome and its functional validation in a Drosophila melanogaster model.

Two Jack-Russell Terrier × Chihuahua mixed-breed littermates with Leigh syndrome were investigated. The dogs presented with progressive ataxia, dystonia, and increased lactate levels. Brain MRI showed characteristic bilateral symmetrical T2 hyperintense lesions, histologically representing encephalomalacia. Muscle histopathology revealed accumulation of mitochondria. Whole genome sequencing identified a missense variant in a gene associated with human Leigh syndrome, NDUFS7:c.535G > A or p.(Val179Met). The genotypes at the variant co-segregated with the phenotype in the investigated litter as expected for a monogenic autosomal recessive mode of inheritance. We investigated the functional consequences of the missense variant in a Drosophila melanogaster model by expressing recombinant wildtype or mutant canine NDUFS7 in a ubiquitous knockdown model of the fly ortholog ND-20. Neither of the investigated overexpression lines completely rescued the lethality upon knockdown of the endogenous ND-20. However, a partial rescue was found upon overexpression of wildtype NDUFS7, where pupal lethality was moved to later developmental stages, which was not seen upon canine mutant overexpression, thus providing additional evidence for the pathogenicity of the identified variant. Our results show the potential of the fruit fly as a model for canine disease allele validation and establish NDUFS7:p.(Val179Met) as causative variant for the investigated canine Leigh syndrome.

PCYT2 deficiency in Saarlooswolfdogs with progressive retinal, central, and peripheral neurodegeneration.

We investigated a syndromic disease comprising blindness and neurodegeneration in 11 Saarlooswolfdogs. Clinical signs involved early adult onset retinal degeneration and adult-onset neurological deficits including gait abnormalities, hind limb weakness, tremors, ataxia, cognitive decline and behavioral changes such as aggression towards the owner. Histopathology in one affected dog demonstrated cataract, retinal degeneration, central and peripheral axonal degeneration, and severe astroglial hypertrophy and hyperplasia in the central nervous system. Pedigrees indicated autosomal recessive inheritance. We mapped the suspected genetic defect to a 15 Mb critical interval by combined linkage and autozygosity analysis. Whole genome sequencing revealed a private homozygous missense variant, PCYT2:c.4A>G, predicted to change the second amino acid of the encoded ethanolamine-phosphate cytidylyltransferase 2, XP_038402224.1:(p.Ile2Val). Genotyping of additional Saarlooswolfdogs confirmed the homozygous genotype in all eleven affected dogs and demonstrated an allele frequency of 9.9% in the population. This experiment also identified three additional homozygous mutant young dogs without overt clinical signs. Subsequent examination of one of these dogs revealed early-stage progressive retinal atrophy (PRA) and expansion of subarachnoid CSF spaces in MRI. Dogs homozygous for the pathogenic variant showed ether lipid accumulation, confirming a functional PCYT2 deficiency. The clinical and metabolic phenotype in affected dogs shows some parallels with human patients, in whom PCYT2 variants lead to a rare form of spastic paraplegia or axonal motor and sensory polyneuropathy. Our results demonstrate that PCYT2:c.4A>G in dogs cause PCYT2 deficiency. This canine model with histopathologically documented retinal, central, and peripheral neurodegeneration further deepens the knowledge of PCYT2 deficiency.

Polioencephalopathy in Eurasier dogs.

BACKGROUND: Polioencephalopathies secondary to inborn errors of metabolism have been described in dogs, but few genetically characterized. OBJECTIVES: Clinically and genetically characterize polioencephalopathy in a family of Eurasier dogs. ANIMALS: Three Eurasier dogs (littermates) presented with early onset movement disorders (9 weeks in 2, 4-6 months in 1). Progressive gait abnormalities were detected in 2 of the dogs, persistent divergent strabismus in 1, whereas consciousness and behavior remained intact in all dogs. One dog was euthanized at 25 months. METHODS: Video footage was assessed in all dogs, and Dogs 1 and 2 had examinations and investigations performed. Whole genome sequencing of Dog 1 and further genetic analyses in the family were performed. A cohort of 115 Eurasier controls was genotyped for specific variants. RESULTS: Episodes were characterized by generalized ataxia, as well as a hypermetric thoracic limb gait, dystonia, and irregular flexion and extension movements of the thoracic limbs. Magnetic resonance imaging of the brain in Dogs 1 and 2 identified symmetrical, bilateral T2 and fluid attenuated inversion recovery hyperintense, T1 hypo to isointense, nonenhancing lesions of the caudate nucleus, lateral and medial geniculate nuclei, thalamus, hippocampus, rostral colliculus and mild generalized brain atrophy. Genetic analyses identified a homozygous mitochondrial trans-2-enoyl-CoA reductase (MECR) missense variant in all 3 dogs, and a homozygous autophagy-related gene 4D (ATG4D) missense variant in Dogs 1 and 2. CONCLUSIONS AND CLINICAL IMPORTANCE: We describe a presumed hereditary and progressive polioencephalopathy in a family of Eurasier dogs. Further research is needed to establish the role of the MECR gene in dogs and the pathogenic effects of the detected variants.

Single-cell profiling of bronchoalveolar cells reveals a Th17 signature in neutrophilic severe equine asthma.

Severe equine asthma (SEA) is a complex respiratory condition characterized by chronic airway inflammation. It shares many clinical and pathological features with human neutrophilic asthma, making it a valuable model for studying this condition. However, the immune mechanisms driving SEA have remained elusive. Although SEA has been primarily associated with a Th2 response, there have also been reports of Th1, Th17, or mixed-mediated responses. To uncover the elusive immune mechanisms driving SEA, we performed single-cell mRNA sequencing (scRNA-seq) on cryopreserved bronchoalveolar cells from 11 Warmblood horses, 5 controls and 6 with SEA. We identified six major cell types, including B cells, T cells, monocytes-macrophages, dendritic cells, neutrophils, and mast cells. All cell types exhibited significant heterogeneity, with previously identified and novel cell subtypes. Notably, we observed monocyte-lymphocyte complexes and detected a robust Th17 signature in SEA, with CXCL13 upregulation in intermediate monocytes. Asthmatic horses exhibited expansion of the B-cell population, Th17 polarization of the T-cell populations, and dysregulation of genes associated with T-cell function. Neutrophils demonstrated enhanced migratory capacity and heightened aptitude for neutrophil extracellular trap formation. These findings provide compelling evidence for a predominant Th17 immune response in neutrophilic SEA, driven by dysregulation of monocyte and T-cell genes. The dysregulated genes identified through scRNA-seq have potential as biomarkers and therapeutic targets for SEA and provide insights into human neutrophilic asthma.

Naturally occurring canine laminopathy leading to a dilated and fibrosing cardiomyopathy in the Nova Scotia Duck Tolling Retriever.

Dilated cardiomyopathy (DCM) is characterized by decreased systolic function and dilation of one or both ventricles, often leading to heart failure or sudden death. Two 10-month-old sibling Nova Scotia Duck Tolling Retrievers (NSDTR) died acutely with evidence of dilated cardiomyopathy with myocardial fibrosis. Association analysis using two cases and 35 controls identified three candidate regions homozygous in the two cases. Whole genome sequencing identified a frameshift deletion in the LMNA gene (NC_049228.1:g.41688530del, NP_001274080:p.(Asp576ThrfsTer124)). Three retrospectively identified NSDTRs with sudden death before 2 years of age and severe myocardial fibrosis were also homozygous for the deletion. One 5 year old with sudden death and myocardial fibrosis was heterozygous for the deletion. This variant was not identified in 722 dogs of other breeds, nor was it identified to be homozygous in 784 NSDTR. LMNA codes for lamin A/C proteins, which are type V intermediate filaments that provide structural support to the nuclear membrane. In humans, LMNA variants can cause DCM with sudden death as well as diseases of striated muscles, lipodystrophy, neuropathies, and accelerated aging disorders. This frameshift deletion is predicted to affect processing of prelamin A into lamin A. Pedigree analysis in the NSDTR and functional evaluation of heterozygotes is consistent with a predominantly recessive mode of inheritance and possibly low penetrance in heterozygotes in contrast to people, where most pathogenic LMNA variants are dominantly inherited.

Independent COL17A1 Variants in Cats with Junctional Epidermolysis Bullosa.

Epidermolysis bullosa (EB), characterized by defective adhesion of the epidermis to the dermis, is a heterogeneous disease with many subtypes in human patients and domestic animals. We investigated two unrelated cats with recurring erosions and ulcers on ear pinnae, oral mucosa, and paw pads that were suggestive of EB. Histopathology confirmed the diagnosis of EB in both cats. Case 1 was severe and had to be euthanized at 5 months of age. Case 2 had a milder course and was alive at 11 years of age at the time of writing. Whole genome sequencing of both affected cats revealed independent homozygous variants in COL17A1 encoding the collagen type XVII alpha 1 chain. Loss of function variants in COL17A1 lead to junctional epidermolysis bullosa (JEB) in human patients. The identified splice site variant in case 1, c.3019+1del, was predicted to lead to a complete deficiency in collagen type XVII. Case 2 had a splice region variant, c.769+5G>A. Assessment of the functional impact of this variant on the transcript level demonstrated partial aberrant splicing with residual expression of wildtype transcript. Thus, the molecular analyses provided a plausible explanation of the difference in clinical severity between the two cases and allowed the refinement of the diagnosis in the affected cats to JEB. This study highlights the complexity of EB in animals and contributes to a better understanding of the genotype-phenotype correlation in COL17A1-related JEB.

Special Issue: "Canine Genetics 2".

Wolves were the first animal species to become domesticated by humans, approximately 30,000-50,000 years ago. Human-directed dog breeding over thousands of generations has generated more than 350 recognized breeds displaying surprisingly different phenotypes with respect to morphology, behavior and disease predispositions. The domestication of wolves and the subsequent breeding of dogs can be viewed as one of humankind's oldest and largest genetic experiments and provides us with unique opportunities for research. Dogs have not only become human's best friend but were also described as geneticists' best friend in a past issue of Science. In recognition of the importance of canine genetics, this Special Issue, entitled "Canine Genetics 2", was compiled. It represents a sequel to the former Special Issue "Canine Genetics", which was published in 2019. During the last 15 years, the canine community has heavily relied on a reference genome derived from the female Boxer Tasha. "Canine Genetics 2" includes an article describing a greatly improved version of this important community resource. This Special Issue further contains several reports related to monogenic or complex inherited diseases in dogs. Finally, important aspects of wild canid research, genetic diversity in different populations and canine morphology were investigated.

Identification of novel genetic risk factors of dilated cardiomyopathy: from canine to human.

BACKGROUND: Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis. METHODS: We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts. RESULTS: Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes. CONCLUSIONS: Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

RALGAPA1 Deletion in Belgian Shepherd Dogs with Cerebellar Ataxia.

Several genetically distinct forms of cerebellar ataxia exist in Belgian shepherd dogs. We investigated a litter in which two puppies developed cerebellar ataxia. The clinical signs stabilized at around six weeks of age, but remained visible into adulthood. Combined linkage and homozygosity mapping delineated a 5.5 Mb critical interval. The comparison of whole-genome sequence data of one affected dog to 929 control genomes revealed a private homozygous ~4.8 kb deletion in the critical interval, Chr8:14,468,376_14,473,136del4761. The deletion comprises exon 35 of the RALGAPA1 gene, XM_038544497.1:c.6080-2893_6944+1003del. It is predicted to introduce a premature stop codon into the transcript, truncating ~23% of the wild-type open reading frame of the encoded Ral GTPase-activating protein catalytic subunit α 1, XP_038400425.1:(p.Val2027Glnfs*7). Genotypes at the deletion showed the expected co-segregation with the phenotype in the family. Genotyping additional ataxic Belgian shepherd dogs revealed three additional homozygous mutant dogs from a single litter, which had been euthanized at five weeks of age due to their severe clinical phenotype. Histopathology revealed cytoplasmic accumulation of granular material within cerebellar Purkinje cells. Genotyping a cohort of almost 900 Belgian shepherd dogs showed the expected genotype-phenotype association and a carrier frequency of 5% in the population. Human patients with loss-of-function variants in RALGAPA1 develop psychomotor disability and early-onset epilepsy. The available clinical and histopathological data, together with current knowledge about RALGAPA1 variants and their functional impact in other species, suggest the RALGAPA1 deletion is the likely causative defect for the observed phenotype in the affected dogs.

SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy.

An 8-month-old female Lagotto Romagnolo dog was presented for a 1-month history of an initial severe reluctance to move, rapidly progressing to a marked stiff gait and progressive muscular weakness and evolving to tetraparesis, which persuaded the owner to request euthanasia. A primary muscle pathology was supported by necropsy and histopathological findings. Macroscopically, the muscles were moderately atrophic, except for the diaphragm and the neck muscles, which were markedly thickened. Histologically, all the skeletal muscles examined showed atrophy, hypertrophy, necrosis with calcification of the fibers, and mild fibrosis and inflammation. On immunohistochemistry, all three dystrophin domains and sarcoglycan proteins were absent. On Western blot analysis, no band was present for delta sarcoglycan. We sequenced the genome of the affected dog and compared the data to more than 900 control genomes of different dog breeds. Genetic analysis revealed a homozygous private protein-changing variant in the SGCD gene encoding delta- sarcoglycan in the affected dog. The variant was predicted to induce a SGCD:p.(Leu242Pro) change in the protein. In silico tools predicted the change to be deleterious. Other 770 Lagotto Romagnolo dogs were genotyped for the variant and all found to be homozygous wild type. Based on current knowledge of gene function in other mammalian species, including humans, hamsters, and dogs, we propose the SGCD missense variant as the causative variant of the observed form of muscular dystrophy in the index case. The absence of the variant allele in the Lagotto Romagnolo breeding population indicates a rare allele that has appeared recently.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Identification of an ADAMTS2 frameshift variant in a cat family with Ehlers-Danlos syndrome.

We investigated 4 European domestic shorthair kittens with skin lesions consistent with the dermatosparaxis type of the Ehlers-Danlos syndrome, a connective tissue disorder. The kittens were sired by the same tomcat but were born by 3 different mothers. The kittens had easily torn skin resulting in nonhealing skin wounds. Both clinically and histologically, the skin showed thin epidermis in addition to inflammatory changes. Changes in collagen fibers were visible in electron micrographs. The complete genome of an affected kitten was sequenced. A one base pair duplication leading to a frameshift in the candidate gene ADAMTS2 was identified, p.(Ser235fs*3). All 4 affected cats carried the frameshift duplication in a homozygous state. Genotypes at this variant showed perfect cosegregation with the autosomal recessive Ehlers-Danlos syndrome phenotype in the available family. The mutant allele did not occur in 48 unrelated control cats. ADAMTS2 loss-of-function variants cause autosomal recessive forms of Ehlers-Danlos syndrome in humans, mice, dogs, cattle, and sheep. The available evidence from our investigation together with the functional knowledge on ADAMTS2 in other species allows to classify the identified ADAMTS2 variant as pathogenic and most likely causative variant for the observed Ehlers-Danlos syndrome.

Case report: Sacral agenesis in two boxer dogs: clinical presentation, diagnostic investigations, and outcome.

Two boxer dogs from the same litter were presented at 3 months of age for urinary and fecal incontinence. Both dogs had an abnormal tail consisting of a small stump, an atonic anal sphincter, and absent perineal reflex and sensation. Neurological evaluation was indicative of a lesion of the cauda equina or sacral spinal cord. Radiology and CT scan of the spine displayed similar findings in the two dogs that were indicative of sacral agenesis. Indeed, they had 6 lumbar vertebrae followed by a lumbosacral transitional vertebra, lacking a complete spinous process, and a hypoplastic vertebra carrying 2 hypoplastic sacral transverse processes as the only remnant of the sacral bone. Caudal vertebrae were absent in one of the dogs. On MRI, one dog had a dural sac occupying the entire spinal canal and ending in a subfascial fat structure. In the other dog, the dural sac finished in an extracanalar, subfascial, well-defined cystic structure, communicating with the subarachnoid space, and consistent with a meningocele. Sacral agenesis-that is the partial or complete absence of the sacral bones-is a neural tube defect occasionally reported in humans with spina bifida occulta. Sacral agenesis has been described in human and veterinary medicine in association with conditions such as caudal regression syndrome, perosomus elumbis, and Currarino syndrome. These neural tube defects are caused by genetic and/or environmental factors. Despite thorough genetic investigation, no candidate variants in genes with known functional impact on bone development or sacral development could be found in the affected dogs. To the best of the authors' knowledge, this is the first report describing similar sacral agenesis in two related boxer dogs.

PKD1 Nonsense Variant in a Lagotto Romagnolo Family with Polycystic Kidney Disease.

A female Lagotto Romagnolo dog with polycystic kidney disease (PKD) and her progeny, including PKD-affected offspring, were studied. All affected dogs appeared clinically inconspicuous, while sonography revealed the presence of renal cysts. The PKD-affected index female was used for breeding and produced two litters with six affected offspring of both sexes and seven unaffected offspring. The pedigrees suggested an autosomal dominant mode of inheritance of the trait. A trio whole genome sequencing analysis of the index female and her unaffected parents identified a de novo heterozygous nonsense variant in the coding region of the PKD1 gene. This variant, NM_001006650.1:c.7195G>T, is predicted to truncate 44% of the open reading frame of the wild-type PKD1 protein, NP_001006651.1:p.(Glu2399*). The finding of a de novo variant in an excellent functional candidate gene strongly suggests that the PKD1 nonsense variant caused the observed phenotype in the affected dogs. Perfect co-segregation of the mutant allele with the PKD phenotype in two litters supports the hypothesized causality. To the best of our knowledge, this is the second description of a PKD1-related canine form of autosomal dominant PKD that may serve as an animal model for similar hepatorenal fibrocystic disorders in humans.

Autosomal recessive hyposegmentation of granulocytes in Australian Shepherd Dogs indicates a role for LMBR1L in myeloid leukocytes.

Pelger-Huët anomaly (PHA) in humans is an autosomal dominant hematological phenotype without major clinical consequences. PHA involves a characteristic hyposegmentation of granulocytes (HG). Human PHA is caused by heterozygous loss of function variants in the LBR gene encoding lamin receptor B. Bi-allelic variants and complete deficiency of LBR cause the much more severe Greenberg skeletal dysplasia which is lethal in utero and characterized by massive skeletal malformation and gross fetal hydrops. HG phenotypes have also been described in domestic animals and homology to human PHA has been claimed in the literature. We studied a litter of Australian Shepherd Dogs with four stillborn puppies in which both parents had an HG phenotype. Linkage analysis excluded LBR as responsible gene for the stillborn puppies. We then investigated the HG phenotype in Australian Shepherd Dogs independently of the prenatal lethality. Genome-wide association mapped the HG locus to chromosome 27 and established an autosomal recessive mode of inheritance. Whole genome sequencing identified a splice site variant in LMBR1L, c.191+1G>A, as most likely causal variant for the HG phenotype. The mutant allele abrogates the expression of the longer X2 isoform but does not affect transcripts encoding the shorter X1 isoform of the LMBR1L protein. The homozygous mutant LMBR1L genotype associated with HG is common in Australian Shepherd Dogs and was found in 39 of 300 genotyped dogs (13%). Our results point to a previously unsuspected function of LMBR1L in the myeloid lineage of leukocytes.

Ancient segmentally duplicated LCORL retrocopies in equids.

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.